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Abstract 


Generic one-parameter billiards are studied both classically and quantally. 
The classical dynamics for the billiards makes a transition from regular to 
fully chaotic motion through intermediary soft chaotic system. The energy 
spectra of the billiards are computed using finite element method which has 
not been applied to the euclidean billiard. True generic quantum chaotic 
transitional behavior and its sensitive dependence on classical dynamics are 
uncovered for the first time. That is, this sensitive dependence of quantum 
spectral measures on classical dynamics is a genuine manifestation of quantum 
chaos. 
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Since the pioneering work of McDonald and Kaufman |l|J^ , the study of statistical mea¬ 
sures of quantum spectra has been one of the major themes for quantum chaos. This study 
uncovered universal behaviors for integrable and most of hard chaotic systems. For tran¬ 
sitional systems, i.e., from integrable to hard chaos through soft chaos senario, however, 
few observations have been made BM- Among those observations, Robnik’s analy¬ 
sis showed using level spacing statistics that there is a gradual continuous transition from 
regular to hard chaos. Honig and Wintgen’s analysis [^], on the other hand, discovered 
that level spacing statistics showed extraordinary behavior in the case of the correspond¬ 
ing classical dynamics with complexity, and suggested that the understanding of it require 
a detailed knowledge of the underlying classical dynamics rather than a knowledge of the 
global classical phase space structure. Grasping a generic transitional quantum behavior is 
very important to understand quantum chaos itself which has not been dehned very clearly. 
These themes motivate this study. 

In this Letter we present the results of complete analysis, both classically and quantally, 
of a generic one-parameter family of billiard systems PJT0[|. As the generic one-parameter 
billiards, we studied Dreitlein’s billiards [Q whose boundary consists of two parallel lines 
with their separation as length and arcs of circles of a radius at both extremes. Fig. 1 shows 
a typical Dreitlein’s billiard. If we vary radii of the circlular arcs the system changes from 
a square billiard to the Bunimovitch stadium. We parametrize the billiards by a parameter 
A as follows. 


A = exp (—V B? — 1)^ (1) 

where R is the radius of circular arcs. As we change A from 0 to 1, Dreitlein’s billiards change 
from a square billiard to the Bunimovitch stadium. At some parameter value A = Ac, the 
billiards become a hard chaotic system, i.e., K-system, through soft chaos. Non-analyticity of 
the boundary for Dreitlein’s billiard makes its classical dynamics relatively simple compared 
to the other one-parameter family of billiards with analytic boundaries. This simplicity 
helps to get the complete analysis of its classical dynamics by studying the Birkhoff surface 
of sections, linear stability of periodic orbits, and imbedded island structures. 

We summarize the results of our analysis of classical dynamics of Dreitlein’s billiards 
in the following. As the parameter A increases from 0 to up, the billiard starts to deform 
from a square billiard. Its Birkhoff surface of section (illustrated in Figure 2) shows lots of 
islands contiguous with each other. The islands of stability, surrounding elliptic hxed points, 
are divided by a very thin connected chaotic region. The absence of separatrix motion due 
to non-analyticity of the boundary makes its structure of the surface of section simple, 
i.e., there is only one major connected chaotic region. As we increase A further, most of 
the islands quickly disappear. The most robust islands surround the corresponding period 
two and period four orbits. Our present numerical analysis shows that the island which 
surrounds period four orbit undergoes period doubling bifurcation at A = 0.167. It is in good 
agreement with A = e“^“^ =0.1663... obtained from the study of monodromy matrix. We 

numerically observed that small isolated regular regions from period doubling bifurcation 
disappear in the chaotic sea at about A = 0.198 The most robust island, which surrounds 
period two orbit, is found to disappear at Ac = 0.368. The study of monodromy matrix 
produces Ac = e~^ which gives very excellent agreement. After that, Dreitlein’s billiards 
become fully chaotic system (hard chaotic system) as A increases. 
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Quantum mechanical analysis should solve the Schrodinger equation with Dirichlet 
boundary conditions. This is equivalent to solve the Helmholtz equation, for which £- 
nite element method(FEM) is employed. FEM is, in principle, one of the best method to 
deal with the boundary problems. The advantages of FEM are known to be able to obtain 
both eigenvalues and eigenfuntions at a time, and to provide solutions to many complicated 
problems that would be intractable by other techniques. Our FEM approach is the hrst 
application to the euclidean quantum billiard, to our best knowledge. We normalize the 
area of the billiards to tt so that the mean level spacing is the same for all Dreitlein’s bil¬ 
liards. Since the classical chaoticity comes from the circular arcs of Dreitlein’s billiards, the 
boundary approximation is very important to obtain accurate results. Because of moderate 
curvatures of Dreitlein’s billiards, iso-parametric coordinates give very good approximation. 
According to Henshell’s estimate |^, the typical radial errors are within less than 10“®%. 


By using FEM, we calculate odd-odd parity eigenvalues and eigenfunctions of Dreitlein’s 
billiards up to about 2400th states for each of 46 different parameter values of A. Analysis 
of these data provides the ways to investigate level spacing statistics, spectral rigidities and 
patterns of eigenfunctions. 

In this Letter, we are limited to discuss only the level spacing statistics. The study of 
the other subjects will be followed near future. The hrst 1100 reliable eigenvalues are used 
for each A to obtain accurate level spacing statistics. To study level spacing statistics of 
Dreitlein’s billiards, we tried to £t the data to the Brody distribution which is originally 
suggested by Brody et.al. |]TT| and later used by Robnik 0. The Brody distribution inter¬ 
polates Poisson and Wigner distributions which are characteristic level spacing distributions 
of integral and chaotic spectra, respectively. 


P{s) = as^exp(-6s^+^) 


( 2 ) 


where a and b are determined by normalization conditions, i.e., J^P{s)ds = 1 and 
J^sP{s)ds = 1. The level spacing exponent z/ is a measure of short range interaction 
between energy levels, i.e., level repulsion. One of advantages to use the Brody distribution 
is that we can dehne cumulative distribution ||^. Since it is a smoother function of s than 
F(s), it is more easily htted by data. We htted our data to the cumulative Brody distribu¬ 
tion, hF(s) = Jq P{x)dx, in order to obtain best htted z/. Fig. 3 shows our main results 
of the study. There are three diherent parameter regions which show diherent qualitative 
behaviors. We call these regions Region I, II and III. 

In the soft chaotic region (Region I), i.e., 0 < A< Ac, it clearly shows that the distribution 
is moving from a Poisson to a Wigner distribution. However, it is certainly not a gradual 
continuous transition. As soon as A increases from 0, u increases sharply and reach about 
0.6 when A becomes 0.05. And then it tends to saturate until A becomes 0.175. Notice 
that the lowest value of A in that stretch is obtained surprisingly at A = 0.167, where 
period two orbits lose their stability. There is a big jump in z/ at A = 0.175 where there 
is only one major isolated integrable region left in the classical phase space. After that u 
huctuates again around 0.7 for a while and followed by another big jump at A = 0.36 which 
is again very close to Ac = 0.368 where the system becomes completely chaotic, i.e., K- 
system. The observation in Region I implies that quantum chaotic transition in Dreitlein’s 
billiards is very sensitive to the corresponding classical dynamics. This rehects that there is 
an intrinsic coincidence between quantum chaos and classical dynamics. In the hard chaotic 
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region (Regions II and III), i.e., Ac < A < 1, we can see two different behaviors depending 
on A. For the region Ac <A< 0.85, as we expected, the level spacing distribution is close to 
the Wigner distribution, which is believed to be a universal characteristic behavior of hard 
chaotic system. In Region II, v is bigger than 0.85 except at A = 0.7. This implies strong 
level repulsions between levels. However for the region 0.85 <A< 1, the exponent v shows 
a clear tendency to decrease. This is an unexpected results which needs an explanation. 
The explanation of this behavior will be the main topic of our future work. Finally, we 
would like to mention that the exponent at A = 1 is 0.717 which is in good agreement with 
McDonald and Kaufman’s 0.71 for the Bunimovitch stadium. Through our complete analysis 
of the corresponding classical system, we have found that quantum spectral measures are 
sensitively related to the topological changes of phase space manifolds of the classical system. 

In conclusion in the present Letter, the major hnding in this study is that the quantum 
chaotic transition is not a smooth gradual transition. There are big changes in the distribu¬ 
tion of eigenvalues for small changes in A. It should be noticed that the parameter ranges, 
0.05 < A < 0.167 and 0.175 < A < 0.36, correspond to parameter ranges in which classical 
phase spaces have two major isolated islands and one major isolated island, respectively. 
And the transition values, A = 0.175 and A = 0.36, are very close to the values at which 
period four and period two orbits lose their stability. This means the detailed dynamical 
behavior of a classical system is very important to understand quantum chaotic transition. 
We claim that this is the generic behavior of quantum transition to chaos with soft chaos 
senario. We think this sensitive dependence of quantum spectral measures on classical dy¬ 
namics is a genuine manifestation of quantum chaos. This also fortihes that spectral analysis 
is the legitimate method to study quantum chaos. 

The authors would like to thank the other GTP members, C.S. Park, D.H. Yoon, S.K. 
Yoo, and D.K. Park, for informative discussions. 
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FIGURES 


FIG. 1. A typical shape of the one-parameter billiards. R is the radius of the arc boundaries. 

FIG. 2. The Birkhoff surface of section for A = 10“^'^ after 100,000 iterations. The initial 
point is given at the neighborhood of the boundary of the period two island. 

FIG. 3. Plot for the Brody’s exponent u versus the parameter A 
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